Single-File Diffusion of Atomic and Colloidal Systems: Asymptotic Laws 
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In this work we present a general derivation of the non-Fickian behavior for the self-diffusion of 
identically interacting particle systems with excluded mutual passage. We show that the conditional 
probability distribution of finding a particle at position xt after time t, when the particle was located 
at xo at t = 0, follows a Gaussian distribution in the long-time limit, with variance 2W(t) ~ t 1/2 for 
overdamped systems and with variance 2W(t) ~ t for classical systems. The asymptotic behavior 
of the mean-squared displacement, W(t), is shown to be independent of the nature of interactions 
for homogeneous systems in the fluid state. Moreover, the long-time behavior of self-diffusion is 
determined by short-time and large scale collective density fluctuations. 



Single-file diffusion refers to the motion of particles in 
narrow pores so that the individual particles are unable 
to pass each other and the sequence of particle labels 
remains the same over time. As the mutual passage is 
excluded, the motion of individual particles requires the 
collective motion of many other particles in the same 
direction. This restriction leads to an anomalous behav- 
ior of the self-diffusion for overdamped systems in the 
long-time limit, which has been subject to long-standing 
theoretical investigations |], ||, ||]. For classical system 
normal diffusive behavior has been found for the self- 
diffusion function [Q, ||. So far, rigorous results have 
been derived for the long-time behavior of the conditional 
probability density function (PDF), P(yt, t\yo, 0), of find- 
ing a particle at position Xt at time t, being located at Xq 
at time t = 0, only for one-dimensional hard-rod systems 
using various analytical methods |Q 0, 0). For these 
systems it has been shown that P(yt, t\yo, 0) is Gaussian 
with variance 2W(t) for asymptotically large times. The 
function W(t) is the mean-squared displacement (MSD) 
given by W(t) — A a i 1 / 2 and W{t) — A c t for overdamped 
and classical systems, respectively. The constants A and 
A c depend on thermodynamic properties and short-time 
transport coefficients of the particular systems. Although 
these results have been derived for hard-rod systems only, 
it is generally believed that in systems with prohibited ex- 
change of particle labels the MSD scales like W{t) ~ i 1 / 2 
in overdamped systems and like W(t) ~ t in classical 
systems. 

Experimental evidence confirming the anomalous dif- 
fusive behavior has been found only recently in atomic 
H and colloidal systems The experimental results 
support the conjectured scaling behavior of the MSD 
as given above §. The single filing condition is hard 
to fulfill experimentally for atomic systems, whereas for 
micrometer-sized colloidal particles narrow channels can 
be easily generated || . While atoms follow ballistic mo- 
tion, colloidal particles follow overdamped motion due 
to the presence of the solvent. Moreover, in colloidal 
systems there exist additional correlations due to sol- 
vent mediated hydrodynamic interactions (HI) |l(| . Nev- 
ertheless the experiments, using magnetic colloids with 



long range dipole-dipole pair interactions, show clearly 
the same long-time scaling of the MSD as predicted for 
hard-rod systems with HI neglected. 

So far a general theory for the self-diffusion in homo- 
geneous single-filing atomic and Brownian systems was 
still lacking and will be given in this work. We show that 
the general conjectured scaling of the MSD for single 
filing systems jTj is indeed true and derive asymptotic 
laws for the single-file diffusion of colloidal and atomic 
particles systems with arbitrary interaction potential, 
under the condition that the correlation length between 
the particles is of finite range and the particles interact 
identically. We prove that in limit of times, much larger 
than the time needed for a particle to diffuse the mean 
particle distance, the conditional PDF P(y t ,t\yo,0) is 
normal distributed with variance 2W(t). Moreover, 
we include HI in the case of colloidal systems and give 
explicit expressions for the pre-factor A in the expression 
for the MSD for both classical and overdamped systems. 

Overdamped systems. We consider a one-dimensional 
homogeneous system of N identically interacting Brown- 
ian particles in equilibrium. The system size is of length 
L and assumed to be macroscopically large. The total 
number density of the particles in the system is denoted 
by p — N/L and the mean particle distance is given by 
a = p^ 1 . The interaction potential is chosen to be infi- 
nite at overlap of any two particle centers. Consequently, 
the particles are unable to pass each other. In the over- 
damped limit the time evolution of the particle trajecto- 
ries follow the stochastic equations [fllj 



%(*) 
dt 



= n B [y N (t)]y l (t)+r H (t) 



(1) 



with CIb the so-called backward Smoluchowski operator 
jfij H given by 
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Here, yi(t) is the position of particle i at time t, Dij{y N ) 
is the diffusivity tensor, accounting for solvent mediated 
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hydrodynamic interactions, and U(y N ) is the interaction 
potential among the particles in the system, given in 
units of UbT. We also make use of the supervector nota- 
tion y N {t) = {yi(t), J/jv(i)}- The random noise t]i(t) 
is normal distributed with mean zero and the second mo- 
ment is given by (Vi(t)Vj{t')) = 2D i3 (y N )6(t - t') . We 
emphasize that the stochastic equations have to be inter- 
preted in the Ito sense fQjfl . The basic quantity we want 
to calculate is the MSD, W(t), given by 



w{t) = \(Mt)- yi m 2 ) 



(3) 



for asymptotically large times t. The index i is arbi- 
trary but chosen such that the particle under considera- 
tion is located sufficiently far from the walls for the time- 
interval considered. The brackets, ( . }, indicate an aver- 
age with respect to the joint probability density function, 
P(y^ , t\yg , 0)P eq (yo), of finding the system at position 
y^ at time t and at position yQ at time t — in con- 
figuration space. Here, P eq {y N ) = exp[—U(y N )]/ZN is 
the Boltzmann distribution and Zn the configurational 
integral. The conditional PDF, P{y^ , t\yQ , 0), satisfies 
the differential equation 



atP(y?,t\y»,0) = n\y?]P(y?,t\y",0) 



with initial condition 



limP(y?,t\y«,0) 



(4) 



(5) 



In Eq. (Q), the operator f2 is the Hermitian conjugate of 
the operator Qb ■ As we assume the system to be ergodic 
and in equilibrium, the conditional PDF P(y^ ,t\y^ ,0) 
will approach the Boltzmann distribution, P eq (y N ), in 
the long-time limit. Hence, lim^oc W(t) is of quadratic 
order in the system size. As W{t) is a monotonically 
increasing function of time, the major contribution to 
the MSD at large times is given by particle trajectories 
whose end to end distances are much larger than the 
mean particle distance, a. In the following we look for 
an expression for the MSD in the limit t >r = a 2 /Ds, 
with Ds the short-time self-diffusion coefficient of a free 
diffusing particle. 

To simplify the calculation for W(t) we introduce an 
auxiliary stochastic function x(t), defined by yi(t) < 
x(t) < and the microscopic density p(y,t) — 

Ef 6 (V - % definition we have W(t) = W(t) + 

0{a) with W(t) = l/2([x(t) - x{0)] 2 ). Consequently, 
the long-time behavior of W(t) and W(t) for infinite sys- 
tem size are related by 



lim lim 



W(t) 
W(t) 



1 



(6) 



In order to calculate W(t) we have to derive an explicit 
expression for the conditional PDF, P(xt,t\xo, 0), of find- 
ing trajectory of x{t) at position xt at time t when it was 
located at xq at time t = 0. 



To this end, it is useful to introduce another stochastic 
variable, yi+u{t) < z(t) < yi + M+i(t), and take M so 
large such that \x(t')z(t")) = for all times t',t" < t. 
By homogeneity of the systems, x{t) and z{t) are stochas- 
tically equivalent variables. Now, the conditional PDF, 
P r (xt — z t , t\xo — zo, 0), to find the two trajectories, x(t) 
and z(t), separated a distance x± — z t at time t when they 
were separated a distance xq — zq at time t — is then 
given by 



P r (Ax t - Az t ,t) = J P(Ax t +a,t)P(Az t + a,t)da 

(7) 

Here, Ax t — x t — x$, Az t — z t — zq and we employed 
translationary invariance and stationarity for the condi- 
tional PDF, P(Ax t ,t) = P(xt,t\x ,0). From the knowl- 
edge of P r (Axt — Azt, t), one can determine uniquely the 
conditional PDF P(Ax t ,t) by making use of the convo- 
lution theorem. 

In the following we are going to determine the function 
P r (Axt — Azt,t) from the stochastic properties of the 
microscopic density, p(y,t). To this end, we introduce 
the functional 

h(x u z t \x ,Zo) := p(y,t)dy-j p(y\0)dy', (8) 

together with the general condition 

lim lim Prob ( (x(t) - x(0)) 2 > W(t) 1+e 1=0, 

(9) 

with e > arbitrary small. Consequently, we get by 
definition of x(t), y(t) 



h(x(t),z(t)\x(0),z(0)) = 0. 



(10) 



The conditional PDF, P r (Axt — Az t , t), is then given in 
long-time regime by 



P r (Ax t - Az t ,t) 



lim lim lim 

l— »oo t»r N,L—>oo 



Si [h(xt,z t \x ,z )] 
dh(xt,zt\xo,z ) 



d{x t - z t ) 



(11) 



Here, t^>r indicates the asymptotic limit for large times 
and 5[z] — lim^oo 5i[z] is the delta distribution in the 
Dirac sense. The reason why we take the asymptotic 
limit t 3> r before taking the limit I — > oo is because the 
conditional probability density, P r (Axt — Azt, i), can not 
be normalized in a simple way if / — > oo is taken first. By 
definition of x(t) there corresponds an infinite set of re- 
alizations of p{y,t) to one given realization of x(t) which 
leads to divergent contributions for (6[h(xt, zt\xo, zq)])- 
Performing the asymptotic analysis for long times first, 
we expect to loose all informations on scales of the mean 
particle distance and hence there is some chance that the 
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same form of normalization is valid, as used in Eq. (|1] 
for a finite set of realizations of p(y, t) corresponding to 
one realization of x(t). We emphasize that the right- 
hand-sides of Eq. (^) and Eq. ( pi] ) are only stochastical 
equivalent in the long-time limit, as in Eq. ( pi] ) we get 
additional fluctuations because in Eq. (||) there enters 
no information that the amount of trajectories of x(t) 
traveling infinitely far during a finite time is of measure 
zero. This additional fluctuations do not contribute in 
the long-time limit, as can be shown by inserting in Eq. 
( 14) any value between the upper bound, as given in Eq. 
(9|), and the lower bound, x(t) — x(Q) = 0, for x(t). 

The conditional PDF for the relative motion, P r (Axt — 
Axo,t), can be rewritten in integral form 



* exp 



E 



,(«) 



e/(2i) 



(12) 



with the n-th cumulant of h(xt, z t \xo, Zo). In the 
following we show that in the limit ( > r only the first 
two cumulants contribute to P r (Aa; t — Azt,t). The first 
cumulant reads simply 



= (h(x t , z t \x ,z )) = p [(x t - z t ) - (x - z )] . (13) 



P r (Axt — Az t , t) = lim lim 



1 



d 



2wi^d(x t - z t ) 



The second cumulant can be expressed as 



,( 2 ) = 



(h(x t ,z t \x ,z ) 2 ) - {h{x t ,zt\xo,za)Y 



j -^^2S(q,0) - (exp[-iq(x t - x )] + cxp[iq(z t - z )])S(q,t)} dq 



(14) 



Here, we introduced the dynamic structure factor de- 
fined by S(q,t) = {5p(q, t)Sp*(q,0)} /N, with Sp(q,t) = 
J dy exp[i q y] (p(y, t) — p) and made use of translationary 
invariance of the system. By inspection of Eq. ( pl| ) we 
find that the dominating contributions to the integral at 
large times come from values q <C a -1 . 

To evaluate S(q,t) in the limit of small values of g, 
we make use of the Mori-Zwanzig projector operator for- 
malism 1 14, pi. For this purpose we introduce the pro- 
jection operators P = \p(q,Q)) eg {p(q,0)\/(N S(q,0)) and 
Q = 1 — P. Here, we make use of the Dirac notation and 
the brackets, {-) e q, indicate averaging with respect to 
Boltzmann weight. The time derivative of the dynamic 
structure factor can be written within this formalism as 



dS(q,t) 

dt 



-q 2 D e (q)S(q,t) 



(15) 

with an effective wavevector dependent dif- 
fusion function, D e (q) = H(q)/S(q) and 
the so-called memory function M (q, t) = 
(p(q,0)Q B \Q exp[Qn B Qt}Q\n B p(q,0))e q @. Here, 
H{q) = (p(qMn B \p(q,0))eq/N is the hydrody- 
namic function |T^]. Detailed calculations show that 
M(q,t) = C(q 4 ) if HI is neglected, or taken into account 
in far-field approximation, Dij(y N ) = P>ij{yi — yj) 
p0| . The full treatment of HI in the memory function 
gives contributions of order q 2 , but the correction to 
S(q,t) are only within a few percent, even for highly 
concentrated suspension with short-range interactions 



p7| . For the colloidal systems used in the experiments 
of Rcf. ||, the far- field approximation of Dij(y N ) has 
been shown to be an excellent one Jl^] . 

For identical Brownian particles and HI neglected or 
treated on pairwise additive level the solution of Eq. jl5| ) 
for S(q, t) is given by 



S{q,t) = S(q,0) exp[-q 2 D e t] + 0{qH) 



(16) 



with D e = linig^o D e (q). Substituting Eq. (|l6|) into 
Eq. (pl|), the second cumulant reads, in the asymptotic 
long-time regime, with S — S^OjO), 



DH 



1/2 



(17) 



Here, the dependence of on x t — Xq is of order 
o(t e ) due to the condition Eq. (||). The higher cumu- 
lants are given by n ( - 2n+1 ^ = and = 0(i x l 2 ) for 
n G {1,2,...}. In order to prove that the odd cumu- 
lants higher than one are exactly zero, we have made use 
of the backward operator to be Hermitian with respect 
to the Boltzmann weighted inner product, (fQ B g) eq = 
(g£l B f) eq , and invariance of the system with respect to 
translation and space inversion. The asymptotic scaling 
of the even cumulants can be found with the help of the 
mean value theorem of integration and the fact that the 
correlation length between particles can be taken finite 
for physical systems in the fluid state. For a 2n-point 
correlation function one uses the mean- value theorem for 
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2n — 2 integrals and then proceeds along the same line as 
for the second cumulant. 

Using the explicit expressions for the cumulants we find 
fromEq. @ 



Pr(Ax t - A**,*) 



lim — 

I — >oo 27T 



d£ exp 



(P/(2 5)) 1/2 
2( 7 rL> e t) 1 /4 



exp 



-£7(20 

p(Aa; t - Az t ) 2 



\S(D e t/w) 1 / 2 

(18) 



From Eq. (^) we can derive finally an expression for the 
conditional PDF for x(t), 



P(x t ,t\x ,0) 



iV2 



2 5 1 /2( 7rj Det)l/4 



exp 



/o(x t - £o) 2 



4S , (£>%/ 7 r) 1 /2 
(19) 

Hence, the mean-squared displacement in the asymptotic 
limit is given by 



W(t) = -{DH/n) 1 ' 2 + o{f' 2 ) 
P 



for t > t (20) 



The effective diffusion constant D e can be experimentally 
determined by a short-time measurement of S(q, t) at 
q -c a -1 Jh|. For times, short enough that each individ- 
ual particle does not 'feel' the presence of the other parti- 
cles by direct interactions, the integral over the memory 
function in Eq. ( |l5| ) does not contribute and thus there 
exists a unique relationship between S(q, t) and D e as 
given in Eq. (|h|) with the second term on the right-hand- 
side neglected. It is interesting to see that a long-time 
self-diffusion property like the MSQ can be determined 
to high accuracy from the short-time collective behavior 
of the system. This fact is a direct consequence of the 
single-filing condition imposed on our system. 
Classical Systems. So far we have considered only inter- 
acting Brownian particles in the overdamped limit. For 
classical particle systems following Liouville dynamics, 
like atoms or molecules in the high temperature limit, 
we can essentially follow the same route as for the over- 
damped systems with the Liouville operator substituted 
for the Smoluchowsky operator. The analysis follows the 
same line as for overdamped systems with the result that 
the conditional PDF for classical systems is again Gaus- 
sian to leading order in time but normal diffusive, 



P(x t ,t\x ,0) = (4^^ s t)- 1 / 2 exp 



(x t - x ) 2 



4D s t 

(21) 

with self-diffusion constant D s = v (1 — dp)/(2p) 
for hard-rod systems (c.f. Rcf(|4|, ||)) and D s = 
c s S Cv I (2 Cp p)^ 1 for classical systems in the hydrody- 
namic regime |20|. Here, v is the average velocity of the 



particles, d is the rod length, c s is the speed of sound 
and Cy, Cp are the specific heats at constant volume 
and pressure, respectively. If the low-frequency sound 
modes are suppressed, like in the presence of a randomiz- 
ing background (e.g. porous media), we expect the MSD 
to scale like W(t) ~ t 1 / 2 ||. This interesting property 
of single-filing systems will be subject of a forthcoming 
paper. In conclusion we have derived a general theory 
for the asymptotic behavior of the MSD in single-filing 
systems. We have shown that W(t) is determined by 
the short-time and large-scale collective behavior of the 
particles in the system. This relation holds true for any 
kind of interaction potential as long as correlation length 
between the particles is finite, which is true for any phys- 
ical system in the fluid state. Moreover, this relation is 
unique in the case of classical systems and in the case 
of overdamped systems with HI taken into account in 
far-field approximation. 

We thank C. Bechinger, J.K.G Dhont and G. Nagele 
for helpful discussions. 
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